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The aim of this paper is to show that within the open-system framework the sum-over-modes 
approach a la Casimir [T| leads to the Lifshitz formula for the Casimir free energy. A general 
result applicable to arbitrary geometries is obtained through the use of Ford, Lewis, & O'ConnelPs 
remarkable formula [5J [3] . Additionally, we address the possibility for obtaining the Casimir energy 
as a sum over complex "modes" . We show in this case that the standard sum-over-modes formula 
must be suitably generalized to avert unphysical complex energies. Finally, we apply our results to 
several standard examples. 

PACS numbers: 



I. INTRODUCTION 

In 1948 H. Casimir predicted that two conducting par- 
allel planes placed in vacuum attract each other. In 
Casimir's derivation this attractive force originated from 
the change in the zero-point energy of the quantum elec- 
tromagnetic field due to the presence of the conductors. 
Indeed, the boundary conditions placed on the field by 
the plates significantly changes the density of states of 
the vacuum fluctuations resulting in a disequilibrium of 
the vacuum field radiation pressure, and therefore to a 
force. 

In his 1948 calculation, Casimir compared the zero- 
point energy of the electromagnetic modes vibrating in- 
side the cavity formed by the two planes with zero-point 
energy of the free field PQ. Mathematically this is equiv- 
alent to the following expression 



E(L) = 



E 

L K 



1 L 



(1) 



where we have introduced the symbol [• • • ]^ which indi- 
cates the difference between energy for a cavity of size L 
(finite separation between the planes) and the energy of 
a cavity with infinite length (infinite separation between 
the planes) . The function wk is the dispersion relation of 
the n-th cavity mode and the symbol K denotes the col- 
lection of good quantum numbers. For the parallel plane 
geometry, 
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where p is the field polarization (TE and TM in this 
case), k denote the transverse wave vector (components 
parallel to the surface) and A is the area of the planes. 
In pQ the conducting planes were assumed to be perfectly 
reflecting so that w K = w£(k) = <V|k| 2 + {nir/L) 2 . The 
generalization to planes made of real materials was un- 
dertaken by E. Lifshitz [4] in 1955. In this seminal pa- 
per Lifshitz adopted a different approach from Casimir 



bearing a strong similarity to London's derivation of the 
van der Waals force [5]. In distinction to Casimir, Lif- 
shitz focused on the electromagnetic (EM) field gener- 
ated by the fluctuating electric polarization within the 
two planes. On the surface this "matter-centric" point of 
view seems drastically different from the "field-centric" 
approach taken by Casimir, although the former is able 
to reproduce the latter in some special limits (conducting 
plates at sufficiently large distance). The difference be- 
tween these approaches leads to the question about their 
connection. 

The first step toward the answer was taken by van 
Kampen in 1968 [3] who showed that the near field limit 
of the Lifshitz formula can be expressed as a sum over the 
zero-point energies of the cavity's surface modes. How- 
ever, the first complete mathematical proof of the con- 
nection between Lifshitz and Casimir's approaches must 
be attributed to Schram [7] who showed the equivalence 
of the fully retarded Lifshitz formula and a sum over cou- 
pled cavity modes. The main limitation of Schram's work 
is that the medium composing the planes is allowed to 
be dispersive in a very general way but not dissipative 
so that the mode frequencies of the EM field are real 
functions. Such a limitation is deeply connected with 
the inadequacy of the expression ([I]) when dissipation is 
present. The dispersion relation, wk, is a complex func- 
tion of K for a dissipative medium and the naive adoption 
of (JlJ to compute the Casimir energy results in a complex 
quantity. 

Here we show that the open systems paradigm pro- 
vides the missing link to demonstrate the equivalence 
of the sum over modes approach and the Lifshitz for- 
mula leading to the Casimir effect between dissipative 
media. The paper is organized as follows. In Section [IT] 
we first recall the essential ingredients of the open quan- 
tum systems theory and briefly discuss one of the most 
used and discussed models to describe dissipation within 
quantum mechanics. Starting from there we generalize 



in Section III Schram's approach to include dissipative 
media. In Section [W] we extend the analysis to general 
geometries, elucidating some of the common misconcep- 
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tions about the connection between Casimir's and Lif- 
shitz's approaches. Finally, in the last section we discuss 
some simple examples and their connection to results al- 
ready found in the literature. 



II. OPEN SYSTEMS FOR THE STUDY OF THE 
CASIMIR EFFECT 

The open systems paradigm can be characterized by 
its subdivision of the total system into two interacting 
components. The first is the system, representing the de- 
grees of freedom we wish to follow, and the second is the 
environment, which is composed of all of the remaining 
parties in the microscopic theory. The total system is 
closed and assumed to be in thermal equilibrium. How- 
ever, if our primary interest is in the system we can for- 
mally integrate out the environment degrees of freedom 
rendering the system open. At this level the environment 
variables make no explicit appearance in the reduced dy- 
namics of the system, and although "hidden" the envi- 
ronment can exchange energy with the system. For the 
treatment of the Casimir effect in this paper, the EM 
field will be treated as the system, and the environment 
is composed of two parties, non-dissipative matter which 
characterizes the local electric polarization within a di- 
electric body, and a bath whose purpose is to introduce 
dissipation into the remaining degrees of freedom. 

Recently, motivated by a longstanding controversy on 
the role of dissipation in the Casimir effect (for a review of 
the controversy see [U [§] and for the related experiments 
see [H2ITT]), several authors have derived the Lifshitz for- 
mula [3] within the open systems approach [T2Tll7| . The 
reason that these different approaches lead to the same 
result can be understood by recalling that the starting 
point for Lifshitz 's calculation is one of the main results 
of Rytov's theory of stochastic electrodynamics [T5]. If 
j = — ituP is a zero average fluctuating current (P is the 
polarization field) the second equation of Lifshitz 's paper 
can be written as 
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is the energy per mode at temperature T, and where 
e(w) is the dielectric function. This relation which might 
seem mysterious at first sight can be made clear within 
the framework of linear response theory. Indeed, Eq. |3]) 
follows directly from the fluctuation-dissipation theorem 
(FDT) [I9H2T] for the fluctuating current. 

In thermal equilibrium the FDT describes the station- 
ary exchange of energy between a system and its environ- 
ment through dissipation and fluctuation channels |22j . 



It is the generalization of the Nyquist theorem [53] 
to the quantum-realm and for arbitrary dielectric bodies 
(Nyquist should be credited with anticipating the quan- 
tum form of his relation in his 1928 paper, see Eq. 8 
of [23 ). To give an intuitive feeling for the physics of 
the FDT we will discuss the latter theorem as a specific 
example that, for a resistor in thermal equilibrium with 
its environment, relates the resistor's current fluctuations 
(Johnson's noise) to its resistance and temperature. 

Imagine that a thermal fluctuation induces a small in- 
stantaneous current through a resistor. This energy is 
then dissipated via Ohmic losses and one might guess 
that the temperature of the resistor will be lowered. How- 
ever, I 2 R (Ohmic) heating accompanies this dissipation 
and feeds back into the Johnson's noise. This process 
characterizes the system's approach to equilibrium and 
in equilibrium these competing processes are balanced. 
In this cycle, despite constant and irreversible dissipa- 
tion of thermal-fluctuation-induced currents, no energy 
is lost but only transformed and reabsorbed in the form 
of heat. The quantum version of this phenomenon adds, 
in addition to thermal fluctuations, the zero point fluctu- 
ations, so that the previous cycle can also take place at 
zero temperature. Even at absolute zero energy can be 
exchanged between coupled subsystems [22l I25H27] . Ul- 
timately, any approach leading to ([3]), the open systems 
paradigm historically 28 -32 being the simplest and most 
powerful, will lead to Lifshitz's formula. A careful anal- 
ysis of each step in the derivation [12 16 has, however, 
the great quality of providing an understanding of all of 
the subtleties in the underlying physics. 

Here we take a different approach to the problem 
which, despite being deeply rooted in the system+bath 
paradigm, does not follow Lifshitz's point of view or that 
of recent work. Rather, we will adopt Casimir's approach 
since the modes of the total system will be the central 
object of interest in our calculation. Of course, since we 
will reproduce the Lifshitz formula, both points of view 
must be strongly connected. 

Our approach will be the following; first, we will an- 
alyze the dynamics of the matter + bath environment, 
which we will refer to collectively as the medium, forget- 
ting at first sight that we are interested in dissipation. 
We will do this by solving the equations of motion for 
bath and then the matter which leads to the medium- 
influenced equations of motion for the field. The effects 
of the medium will be accounted for through a permit- 
tivity function for the material and a Langevin forcing 
term. The concert of these two effects describe how the 
field and the medium exchange energy. The nature of this 
energy exchange is roughly described by the Poincare re- 
currence time which quantifies the time required for the 
system to return to a state resembling its initial config- 
uration. For a countable (even infinite) set of coupled 
oscillators, as we are using to describe our environment, 
the Poincare recurrence time is finite, and therefore any 
energy lent to the medium by the field will be returned in 
a finite (possibly very long) time. That means that the 
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flux of energy through the dissipation channel is not uni- 
directional and therefore the medium is non-dissipative, 
in the usual sense, so long as the bath is composed of a 
countable set of coupled oscillators. Our next step will 
be to study the interaction of this macro non-dissipative 
medium with the EM field. At the end of the calcula- 
tion we account for dissipation by allowing the Poincare 
recurrence time to become sufficiently long (infinite) so 
that the dissipation channel becomes unidirectional. In 
this limit, which is implemented by taking the frequen- 
cies of the bath oscillators to form a continuum, energy 
is irreversibly lost to the bath through dissipation. This 
does not imply that the system cannot remain in equi- 
librium. For the case of a resistor the dissipated energy 
heats up the environment which enhances the Johnson's 
noise, these two competing processes eventually balance 
and maintain the system in equilibrium. An analog pro- 
cess will enhance field fluctuations as energy is dissipated 
to maintain thermodynamic equilibrium for the system 
of interest in this paper. 



A. A simple model for a dielectric medium 

To clarify the previous concepts let us consider a simple 
model for a dielectric medium where the optical proper- 
ties are described by a collection of independent harmonic 
oscillators. (The following procedure can be generalized 
for more complicated models, like a medium composed of 
interacting oscillators.). The dynamics of each oscillator 
is described by the equation of motion 



m[x(t) 
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where the frequency ujq describes the elastic response, m 
the mass, and F ext (t) is an external force that we will 
assume to be proportional to the electric field, F ext (t) = 
— eE(f) (e quantifies the coupling to the electric field). 
We will assume that the field can be considered to be 
spatially homogenous over the displacement of the os- 
cillator (long wavelength approximation). The steady- 
state dynamics of the oscillator is accurately described 
by its susceptibility (polarizability) which quantifies the 
(linear) response of the oscillator to an external pertur- 
bation 

- ex(w) = d(oj) = a(u)~E(uj) (6) 
where d(w) is the dipole moment, from which we derive 

e /m 



a{uj) 



LU7{ — UJ' 



(7) 



Given this we can conclude that the dielectric constant 
of the material is given by 
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with uip = Awne 2 /m, where n is the density of oscillators 
per unit of volume. Note that at this stage the dielectric 



function is a real-valued function reflecting the time re- 
versal symmetry of the system (see also the discussion in 
the following sections). 



B. Coupling with the bath 

Now, within this simple model for a dielectric we 
are going to employ the open system approach. Let 
us suppose that each oscillator composing the dielectric 
medium is coupled to a bath of harmonic oscillators. A 
typical Hamiltonian for this system is given by [5^1 130] 
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Note that this is one of many possible Hamiltonians 
which can describe the medium, the main differences are 
due to the specific choice of system-bath coupling. (Other 
examples are given in Refs. [30 , 35] where the bath is pro- 
vided by the EM field and the effect of the electron form 
factor is considered.) However, the relevant features we 
are going to discuss are common to all models. This gives 
equations of motion 
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The formal solution of the second equation is given by 

q» = qf°» + , xM, (13) 



which when combined with the previous equations gives 
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-f(oj) (Langcvin force) 



In the previous equation the prefactor of d(w) is the in- 
verse of the generalized polarizability 
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while the second term on the r.h.s. is a Langevin force 
[2"9"] describing the bath-induced fluctuations of the oscil- 
lator. 
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Some remarks are due before proceeding. The cxpres- is continuous, making an infinite Poincare time possible. 



sion for the polarizability given in Eq.(15) is more gen- 
eral than the specific model considered here and it has 
strong connections with the theory of quantum Langevin 
equations [30 : only the detail of the expression for /j,(uj) 
directly depends on the specific coupling with the bath. 
General physical principles impose that fi(uj) is a posi- 
tive function, which means that (i) it is analytical in the 
upper-half of the complex- frequency plane (causality con- 
dition), (ii) it has a non negative real part at the upper 
boundary of the real axis (to preserve the second law of 
thermodynamics), and (iii) it must satisfy the crossing 
relation = p*(—(*) (reality of the function in the 
time domain). (For further details see for example Ref. 

Note that the generalized polarizability is at this stage 
real and quadratic in lj. However, if we take the limit 
that the bath oscillator's frequencies form a continuum 
we find 



d// 
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where z0 + has been added to the frequency in order to 
enforce causality and we have introduced the bath den- 
sity of modes p Uc (v), with uj c a cut-off frequency such 
as pu c { v > w c) ~ 0- We will assume that the product 
of m(v) and bath density of modes is even in v allowing 
the ^-integral to be extended from — oo to +oo . As a 
particular case, if we assume that the only poles in the 
previous expression are at v = ±(cu + i0 + ) and that the 
integrand vanishes for large imaginary v the integral can 
be evaluated to give 



l m 
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By further choosing the spectral density such that /i(w) = 
7 the generalized polarizability reduces to the well 
known Drude-Lorentz form for the polarizability 



Eq.(15) is still a real quantity with the following proper- 
ties: 

a{ui) = a(— uj) reversibility, (22) 
a(ui) = 0. (23) 

The solutions to l/a(uj) = describe free oscillations 
that can exist in the absence of external forces, i.e. the 
modes of the coupled system (this can be easily seen in 
the case of Eq. ffi where the pole is given by the oscilla- 
tor frequency lJq) [5] . These modes are collective oscilla- 
tions of the whole medium. The tensor product of energy 
eigenstates of the uncoupled system are not eigenstates 
of the total system because, at equilibrium, their degrees 
of freedom are entangled |5H [33] . This means that de- 
spite the fact that each collective mode has a well defined 
fixed energy, this energy will flow back and forth between 
the matter and bath degrees of freedom. These proper- 
ties are general results of linear response theory and also 
apply to more complicated systems where the role of the 
polarizability is played by the the Green's function (see 
next sections). 



In the following we will forget about Eqs.(19) and (20) 



and we will work exclusively with Eq.( 15 1 and the corre- 
sponding real-valued dielectric function 
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where the subscript "iV" indicates that the heat bath is 
made by a countable number of oscillators. This function 
can be directly used in Schram's approach to derive the 
Casimir energy from which the Lifshitz formula can be 
recovered by taking the continuum limit {N disappears) 
for the spectrum of bath oscillators in the final formula. 



III. CASIMIR ENERGY A LA SCHRAM 
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which, when combined with Eq.Q, gives the usual 
Drude-Lorentz model for the dielectric function 



e(w) = 1 - 



;(w + ry) — uj 2 



(20) 



In the time domain and in the absence of an external 
electric field, we have 



d(t)+7d(t) + wgd(t) = — f(t) 



(21) 



which shows how the dipole undergoes dissipation and 
how the heat bath drives the dipole through a Langevin 
force One has to remember, however, that this result 
can only be obtained when the bath oscillator spectrum 



One can easily verify that the dielectric function in 
Eq. (24) is an even function of frequency, i.e. e./v(w) = 
ejv(— uj), which results from the time-reversibility of a 
closed coupled system composed of a discrete number of 
harmonic oscillators. This condition is sufficient to de- 
rive the Casimir energy using a sum-over-mode approach 
similar to the one used by Casimir in his seminal pa- 
per. To do this we will follow the procedure proposed 
by Schram in [7] which avoids problems with branch cuts 
(see Fig[T]). For the development described here it will 
suffice to point out the main steps of [7]; a more thor- 
ough reading of Schram's work [7J is recommended for 
the reader in search of more details. 

Following Schram's proposal our cavity is formed by 
two identical mirrors of thickness d and composed of 
a material whose optical properties are determined by 
ejv(w). The external surfaces are chosen to be perfectly 
reflecting (we put the dielectric mirrors inside a perfectly 



5 



(a) 



4. 



L 



non-dissipative 
limit 



b 


















untab 










i 


o 
o 
d 
Z3 




d 


dissipative limit 





















(c) 



.©■ma 



Figure 1: (Color online) Schematic representation for the im- 
plementation of dissipation in the sum-over-mode approach 
to derive the Casimir effect. Subfigure (a) depicts a cavity 
formed by parallel planes where the optical properties of the 
mirrors are determined from the dynamics of the electric po- 
larization field coupled to a bath of harmonic oscillators. The 
system as a whole is closed and at equilibrium. To avoid 
branch cuts we have introduced two perfectly reflecting mir- 
rors on the outer edges [7J (see text for further explanation). 
As long as the number of bath oscillators is countable the 
modes of the total system are real and the dielectric medium 
does not exhibit dissipation. Subfigures (b) and (c) show that 
the local electric polarization within the dielectric body cou- 
ples to the bath's "hidden" degrees of freedom allowing the 
bath and polarization field to exchange energy. In the non- 
dissipative limit depicted in (b) the polarization field couples 
to a discrete countable number of oscillators. For this case 
energy lost to the bath from the matter will return in a fi- 
nite time. Dissipation is introduced when the number of bath 
oscillators becomes uncountable. This is illustrated in subfig- 
ure (c) where the discrete oscillators shown in (b) are smeared 
into a continuum. In this case dissipation manifests as an ir- 
reversible transfer of energy from the polarization to the bath 
resulting ultimately in the thermalization of the matter. 



reflecting cavity). By imposing the boundary conditions, 
we find that eigenfrequencies of the EM field vibrating 
in this multi-layered planar cavity |36j are given by the 
zeros of the following expressions: 



Djf E (uj,L) — e KL (k + n m taiih[K m d]) 2 
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(herein we measure frequency as a wavevector, effec- 
tively changing variables so that ujjc — > u>) and where 
we choose the branch of the square root to be defined by 
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Figure 2: (Color online) Two different but equivalent choices 
of paths for the calculation of the Casimir (free) energy 
as a sum over modes. The red dots along the imaginary- 
frequency axis indicate the position of the Matsubara's fre- 
quencies (poles of the hyperbolic cotangent in Eq.Q) while 
the black dots along the real-frequency axis are the frequency 
modes of the total field+(matter+bath) system. 



Re k, Re K m > and Im k, Im n rn < 0. The presence of 
perfectly reflecting surfaces on the outermost faces of the 
cavity leads to the appearance of the functions tanh[K m d] 
and coth[K m d] in (25) which render D v N (u),L) even with 



respect to the variable K m . Thus, in a series expansion of 
D p n (lj, L) in powers of K m only even powers of n m will ap- 
pear which guarantees that no square roots involving the 
frequency will appear from the dependence of the func- 
tion D^(lu, L) on K m . Eqs. (25 1 still possess a branch cut 
in the complex frequency plane through dependence on 
the variable k (they are odd in the variable k), but this 
complication can be avoided by multiplying the previous 
functions by k; this extra prefactor does not depend on 
any geometrical parameters of the system and will auto- 
matically be removed in the renormalization procedure 
(see the following). This procedure is, however, neces- 
sary for the consistency of the approach and guarantees 
that Eqs. (25) only depend on powers of the frequency. 



Therefore, the functions kD p n (u, L) do not have branch 
cuts in the complex w-plane, which is a consequence of 
the discretization of the field modes due to the perfectly 
reflecting mirrors on the boundaries [7] and multiplica- 
tion of the dispersion relation by n. 

The zeros uj m (p = TE, TM) of Eqs. ([25]) are the 
modes of the of the coupled system field+ (matter +bath). 
Note that since the dielectric function is real all the 
modes lie on the real frequency axis. The system does 
not dissipate at this stage because the energy can still 
come back after a time corresponding to the Poincare 
time of the system (see Fig[T]). Dissipation will be in- 



6 



eluded in the last step of our derivation by taking the 
Poincare time to infinity. Until this moment we can fol- 
low Schram's derivation and calculate the energy of the 
total system just by adding the energy of each mode us- 
ing the argument of the logarithm theorem |37) 



S N (L) = ^ u(u! K ) 
dz 
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2m 



u(z)d z \n[ K D p N {z,L)}, (27) 



where, u(ui) is the energy per mode defined in Q, and 
C is a path that encloses the right part of the complex- 
frequency plane in the positive sense and displaced a van- 
ishingly small ditance, <5, away from the imaginary axis 
in order to avoid the poles of the hyperbolic cotangent 
(see Fig{2]) . This quantity is divergent because it contains 
by default the energy of the free vacuum plus the energy 
of the modes vibrating inside (or along the surface) of 
each isolated mirror (more on this point in the follow- 
ing) . This energy can be removed (it does not depend on 
the distance between the mirrors) by subtracting from 
the previous expression its value in the asymptotic limit 
L — > oo. Therefore, we define the Casimir energy of our 
system as 



E N (L) = [S N (L)} 
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where it can be seen that the factor k, in the argument 
of the logarithm cancels. The previous function is con- 
vergent as we can see, for example, in the case of p = TE 
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and similarly for p — TM . Now, since the integrand goes 
to zero for \z\ — > oo we can eliminate the part of the 
contour C at infinity and after an integration by parts 
and a change of variables one gets 



E N (L) = V ^ ^u(i£ + 5) In g p N (i£ + 5, L). 



With the help of the following relations 

= -1. = 
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d z u{z) = ~T 2 d T ^, (31) 



zT 



where the volume is fixed when taking the temperature 
derivative, we get the Casimir free energy (from which 
one derives the force) 

F N {L) = k B T Y, ^g p N {it u L). (32) 
p,k,;=o 



In the last step we took the limit <5 — !• using 
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(33) 



then we used the parity of the integrand to eliminate 
the contribution of the principal value. The frequencies 
£l = 2TrkBTl/(hc) are the Matsubara's frequencies and 
the prime in Eq. ( 32 1 indicates the the zeroth term of 



the sum over I has weight 1/2. An equivalent expression 
(which was chronologically obtained first by Lifshitz [4]) 
can be derived using instead the contour Ci encircling 
(a vanishingly small distance, iS±, from the real axis) the 
real-frequency axis (see Fig. [2]) and summing over the 
energy Q or directly free energy per mode 



f(u>) = k B T\n 2sinh 



2k B T 



(34) 



(The previous function has branch points for z = i£i 
and branch cuts in the left part of the complex-frequency 
plane.) The contribution from the small arcs at zero and 
at infinity can be removed when the path wraps tightly 
around the real frequency axis [Si — ¥ 0) and this gives us 
the free energy 



T N {L) = - E / — /(w)I*>4, In G p N (u> + i0 + ,L) 



p.k 



= ^E 





Hcuj 


/ — coth 




Jo 2ir 


2k B T 



Imln g$f(u + iO + ,L). 

(35) 



This expression is, however, of less practical value for ac- 
tual computation as its integrand rapidly oscillates. The 
equivalence of both expressions can be proved by per- 
forming a Wick rotation 0] in the complex frequency 
plane, taking into account the integrand is analytical in 
the upper-half of the plane and that the contribution 
from the integral along the arc at infinity vanishes. 

Before proceeding let us take a step backwards and give 
some physical interpretation of the mathematics. For- 
mally, the interesting part of the renormalization pro- 
cedure is the subtraction from £n(L) its corresponding 
value at L — > oo. To be completely rigorous in the follow- 
ing discussion we should introduce a cut-off function to 
parameterize the divergences. We avoid such formality 
at this stage with the claim that we get the same result 
from a more rigorous derivation. One can show that the 
energy £^{L) for L — ¥ oo can be written as follows (we 
will consider p = TM only but we get something equiv- 
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alent for TE too) 
£ T N M {L) L ^°° LA J 



d 2 k dk z 
(2tt) 2 ~2V 



i(v/k 2 + fc2) (36a) 



do; 



2y I — u(w)Im9 w In (e^(w)fi; + /t m coth[«; m d]) 

(36b) 



We recognize in the line ( 36a I the vacuum energy of the 
free (without mirrors) EM field in the volume of the cav- 
ity (V — LA). In line (36b I if we apply the argument 
principle we get two times the sum of the zero point en- 
ergies corresponding to the solutions of 



epf(uj)K + K m coth[n m d] = 0. 



(37) 



We recognize that these solutions give the plasmons, 
eddy currents, and bulk plasmons plus all other modes 
vibrating inside a single mirror. The factor two in ( 36b ) 



simply means that we are dealing with two identical 
mirrors (everything can be straightforwardly generalized 
to the case of different mirrors). The previous expression 
gives a physical description of how renormalization 
works. 



Going back to (32 1 we are two steps from recovering 



the Lifshitz formula as in his 1955 paper. These two steps 
are two limits: 1) the continuum limit for the spectrum of 
bath oscillators and 2) the bulk approximation. The first 
limit transforms ejv(w) into e(w) and the second limit is 
performed by taking the thickness of the dielectric mir- 
rors, d, to infinity giving 



tanh[« m (i], coth[K m d] d ~ > °°> 1. 



(38) 



Taking these two limits we recover the familiar Fresnel 
reflection coefficients as a function of imaginary frequen- 
cies 

TE/-/-\ K ~ K m , TMl-f\ € (K) K ~ K m 



so that Eq. ( 32 ) reduces to 



e(*£) K + K m' 

(39) 



T(L) = k B T ln[l-r p («e 
p,k,i=o 



2 e -2 KL] 



(40) 



i.e. the Lifshitz formula. In the zero temperature limit 
(J 7 = E) we have 



-2kL~\ 



(41) 



As a last remark of this section let us stress the impor- 
tance of taking the continuum limit for the bath spec- 
trum at the end of the calculation (the following con- 
siderations also apply to the next sections). Before this 



oscillator 
• 



Countable bath modes 



Continuum Countable oscillator+bath modes 

limit 




Figure 3: (Color online) Schematic diagram illustrating how 
dissipation is introduced by taking the continuum limit for 
the spectrum of the bath oscillators. In the system+bath 
paradigm the system (assumed here to be a harmonic oscilla- 
tor for simplicity) is coupled to the bath. The modes of the 
total system (green-red dots) are therefore the hybridization 
of the modes of the system (green dot) and of the bath (red 
dots) and are all described by (countable) real frequencies. At 
this stage it is still possible to analytically continue all related 
expressions from the upper to the lower part of the complex 
frequency plane along a path that does not go through a pole. 
In the continuum limit the mode spectrum becomes dense on 
the real axis making it impossible to access the lower part 
of the complex plane. An analytic continuation, in analogy 
to what happen with branch cuts, leads to an "unphysical" 
Riemann sheet. 



limit is taken all quantities are even in frequency (micro- 
reversibility) and are well-defined along the positive and 
negative imaginary frequency axis. However, this is no 
longer true for the final expression as one can easily verify 
from the Fresnel reflection coefficients. These have, in- 
deed, branch cuts along the negative imaginary frequency 
axis, making an integration along this axis ill-defined. 
This behavior is attributable to the fact that we must 
impose boundary conditions on the bath susceptibility 
to derive a causal permittivity (see the i0 + in Eq. (17)). 
Therefore, 



before taking the dissipative limit, we make 
all necessary manipulations of our integral expression for 
the Casimir energy so that all integrations are in the up- 
per half of the complex plane. 

Also, as was pointed out by Ford, Lewis, and O'Connell 
in [2] and by Sernelius in [38] , the use of the expression for 
the causal permittivity in the lower part of the complex 
plane is inappropriate. Indeed, in the continuum limit 
the modes along the real axis form a cut. Therefore it 
is misleading to use the same expression for the permit- 
tivity (with retarded boundary conditions applied to the 
bath) in the lower half of the complex frequency plane 
(see Fig{3]). An analytic continuation to the lower plane, 
in analogy with what happens with branch cuts, leads 
to an "unphysical" sheet [2J. Within the open system 
approach, the complex poles that appear in this case are 
a mathematical abstraction which signals the manifesta- 
tion of an infinite Poincare time, to which we are used to 
attaching a physical meaning as complex "modes" of the 
system but that in reality have little to do with modes 



since they do not have a well defined energy (see also in 
the following and in Ref. [341139]). 



IV. GENERAL GEOMETRY AND 
APPLICATION OF THE REMARKABLE 
FORMULA 

The previous analysis has showed that it is possible to 
derive the Lifshitz formula using a sum-over-modes ap- 
proach even for a dissipative system. The key point was 
to understand that what we usually call dissipation is 
nothing but the addition of an extra system (bath) with 
an infinite and continuous number of degrees of freedom. 
The addition of a bath modifies the sum-over-modes tech- 
nique by incorporating new modes into our system (noise 
lines) that contribute to the zero point energy of the total 
system. The usual form for the Casimir energy is recov- 
ered by adding the energy of all modes and by taking the 
continuum limit for the spectrum of bath frequencies. 

The previous section ended with the Lifshitz formula 
for the Casimir energy between two parallel planes. In 
this section we show that the same result can be ob- 
tained in a more general framework with the use of the 
Ford, Lewis and O'Connell remarkable formula, follow- 
ing an approach very similar to what was presented in 
[3] . (The connection between dispersion forces and the 
remarkable formula was already recognized by Obcemea 
in [40] .) Here we show that, together with the theory 
of open systems, Ford, Lewis and O'Connell's approach 
allows us to derive the Casimir (free) energy in full gen- 
erality (including dissipation). The final formula can be 
written as follows 



F{L) = -he [ 
Jo 



f°° 




/ — coth 




/o 27T 


2k B T 



Im Tr In 



G(u) 



GooM 

(42) 

where G{uj) is the electromagnetic Green tensor corre- 
sponding to the geometrical configuration (plane-plane, 
sphere-plane, etc). Goo(cj) is the Green tensor of a ref- 
erence configuration with respect to which we measure 
the Casimir energy; as an example, for "outside geome- 
tries" [4*1] . like two parallel planes, this will generally be 
assumed to be the configuration where the two objects 
are infinitely spaced, i.e. where their interaction energy 
vanishes. This definition loses its physical clarity when 
we deal with "inside geometries". In this case it is no 
longer possible to infinitely separate the objects without 
changing the "inside" features of the positioning, mak- 
ing the reference configuration corresponding to Goo(a;) 
more arbitrary. For example, in the case of two concen- 
tric cylinders this Green tensor can correspond to the 
sum of configurations where only one of the cylinders is 
present (sum of the self-energies) which mimics the non- 
interacting configuration in a manner similar to the case 
of two planes placed at infinity. An alternative could be 
the configuration where the two cylinders are not only 
one inside the other but also coaxial. From the energy 



point of view the two "renormalization" procedures lead 
to two different results; in the second case the reference 
energy contains together with the self-energies of the two 
isolated bodies a part of the interaction energy. This 
extra energy, however, does not lead to an extra force 
that, for symmetry reasons must be zero when the two 
cylinders are coaxial (stable or more probably unstable 
equilibrium point). Physically, the difference depends on 
how we understand the Casimir energy: If we define it as 
the energy due to the interaction of two bodies we must 
define it with respect to the self-energies; if it is just the 
work done by the Casimir force we have the second result. 
The difference is only a constant. At zero temperature 
this constant is insignificant but could be more impor- 
tant at finite temperature. We will analyze this point 
with more details in a forthcoming publication and here 
we will limit our considerations to outside geometries. 



A. Equations of motion 

In this section we set the fundamental concepts that 
will lead to Eq. (42). Let us start from the relevant 



macroscopic Maxwell equations written as 

V x V x E(w, r) - w 2 E(w, r) = Attuj 2 P(uj, r) (43) 

where E is the electric field, and P is the matter field 
which quantifies the local electric polarization within the 
interacting objects. To go further we need a model which 
will provide the equations of motion determining the dy- 
namics of the matter field. To this aim let us consider the 
Hamiltonian density for an elastic (hydrodynamic) field 
coupled to a collection of (for simplicity incompressible) 
oscillators (bath) 

This description is just the generalization of the model 
used in the previous section and it can easily be con- 
nected with the well-known model discussed by Huttner- 
Barnett in previous work |42j . The corresponding equa- 
tions of motion are (in Fourier space) 

- ff 2 [/3 2 V 2 + (c 2 -c 2 )]P(c,r) = 

N 

E(w,r)+5^ Wj a {^(r)-P(r)} (44) 



and 



[- W 2 + W 2 ]^(r)=c 2 P(r) 



(45) 



where /? is the compressibility factor and the gj 's are con- 
stants characterizing the coupling strength of each bath 
oscillator to the matter. The compressibility factor intro- 
duces a simple spatially nonlocal component in the elastic 
response of the polarization fluid |43j . This component is 
not included in the dielectric model most often used for 
describing dielectric media but can be easily introduced 
into a Lagrangian approach as the one described in [42] . 
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Solving the second equation for cf>j (r) and inserting its 
solution into (44) we find 



-co 



2v72 



, , 2 9 i , , 2 



P(w,r) 



= E(w,r) +f(w,r) (46) 



where we have defined 



f( W ,r)=£#^tff>,r) 



(47) 



which describes a random fluctuating (Langevin) polar- 
ization. 

The general solution for the previous equation can be 
written as follows 

P(w,r)=P noise (u;,r) + J d 3 r'^ N (ij; r, r') • E(w, r'), 

(48) 

where P no j so (oj,r) is the polarization noise induced by 
the Langevin force f(w,r) 



,(« >r ) = J dV^ (a;; r,r') •%,!•') 



(49) 



and ^jy( w ; r 7 r ') is the Green tensor of Eq.(46). Since 
the polarization and bath field only have support 
within the volume defined by the interacting objects 
X iv( £i 'i r > r ') contains all information about the geom- 
etry of the system. We should point out here that we 
have dropped the homogeneous solution of the polariza- 
tion field. In the dissipative limit the free evolution of the 
polarization is damped in time; at late times, i.e. thermo- 
dynamic equilibrium, the fluctuations of the polarization 
are determined entirely by the bath. 

At this point it is convenient to adopt an operatorial 
(not base dependent) notation and write the solution of 
|i6l as 



Xn [|E) + |f)] 



(50) 



The previous expressions can be recovered by projecting 
the previous equation onto the space-basis, for example 



Viv(w;r,r') = (v\xnW) 
E(w,r) = (r|E). 



(51) 
(52) 



It is implicitly assumed that all quantities depend on the 
same frequency, which will be made explicit if necessary. 
Maxwell's equations (43) together with (50) give the fol- 
lowing equation 



Go" 1 - X n] |E) = |P 



(53) 



where we have defined the inverse free vacuum Green 
tensor as 



(r|GoV) 



<5(r- 



47TW 2 



[-LJ 2 + V x Vx] 



(54) 



Using Eq. ( 53 1 we find 



IE) 



GAr|P r 



Gn ■ Xn\() 



(55) 



where G at (w; r, r') is the Green tensor that contains all 
information about the geometry and the dynamics of the 
coupled EM-matter-bath system. As before, we dropped 
the homogeneous solution as it plays no role at late times 
in the dissipative limit. From the previous relations we 
have that 



Xn 



Gn 



G 



N 



(56) 



or also 



G 



N 



1 — Go ■ xn 
Go ■ 1 — Xn ■ Gq 



Go 
-l 



(57) 



Let us conclude this section with some remarks showing 
the simplicity introduced by this notation. Going a step 
ahead we take the continuum limit for the bath spectrum, 
which makes xn X an d Gn — > G, as well their inverse, 
become complex quantities. From (56), since Gq 1 is still 
a real quantity, we have that 

rr»-ii* _ n-i 
Im x = -Im G- 1 = [ J 2 , . (58) 

One can then show the following relation 

G* ■ (Im x) • G = Im G (59) 

which connects the fluctuation-dissipation theorems of 
the first and second kind HH H5I ■ 



B. Modes and Remarkable formula 

For the next step toward deriving the remarkable for- 
mula (42) we need to establish a direct connection be- 



tween the Green tensor (operator) and the modes of 
the total system. The approach will be a generaliza- 
tion of the discussion presented by Ford, Lewis and 
O'Connell in Ref.[5] who showed that the zeros and poles 



of Dct 



give the modes of the total system and 

the modes of the environment, respectively. In the pre- 
vious expression, the symbol "Det" indicates the deter- 
minant for both the position and the polarization coor- 
dinates. We will keep the discussion general so that the 
results we obtain will give the properties of the suscepti- 
bility in the linear-response theory for arbitrary systems 
within the open systems framework. 

To begin, let us consider the Hamiltonian H for a 
generic total system which is a closed. The Heisenberg 
equation of motion for a generic operator A(r) is given 
by 



d t A(t,r) = -[H,A(<,r)] 



(60) 
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If we define the superoperator H = [H,-]/h, the last 
expression, in Fourier space, can also be rewritten as 



H\A(u)) = -u\A(w)) 



(61) 



where |A(cl>)) is the base-independent form of the oper- 
ator A(cj,r). If Ho is the sum of the free Hamiltonians 
of each subpart of the total system we know that the 
equation 



JY 



N 



3=1 3=1 
N 

^-Y^g^jWr^)) (62) 

3=1 

is identically satisfied if |<^ ree (w)) are given by the bath 
operators undergoing free evolution. Our focus here is 
connected with the fact that we are interested in the 
dynamics of the system driven by the bath. Using a 
Lippmann-Schwinger-like approach it is possible to show 
that |A(w)) must satisfy the following equation 

|A(w)) = [H + oj + i- n r 1 -[Ho+oj + i^-H^lfiuj)) (63) 



where in agreement with Eq. ( 47 1 we have 



|f(")) 



JV 

$J a 3 

3=1 



(64) 



To enforce causality we add the infinitesimal imaginary 
factor ir\ (r\ > 0) to the frequency. At the end of the 
calculation we take 77 —> 0. This superoperatorial rela- 
tion shows that each operator of our system satisfies a 
Langevin-like equation. When \A(u)) — |E(w)) and Hy 
is the Hamiltonian of the total EM ficld-matter-bath sys- 
tem from Eq.(55| we get 



[H t + lo + ir,] ■ [H + uj + ir)} ■ H Q 2 = G N ■ X N- (65) 

Similarly, we can consider the dynamics of the bath- 
influenced matter without coupling to the EM field. 
Thus, we choose |A(w)) = |P(w)) and replace H with the 
Hamiltonian for the coupled matter-bath system, Hm- 
By using Eq.([50]) we find 

[H M +oj + irfT 1 -[H q +oj + ir,] ■ H^ 2 = X n- (66) 

From the previous expressions we can deduce 

G N = [H T + u + i-qf 1 ■ [H M + u> + irj\. (67) 

If the total system is linear its Hamiltonian can always 
be diagonalized in the following form 



H = ^2 f^K (B t (K)B(K 



K>0 



(68) 



where K is the collection of quantum numbers that char- 
acterize its eigenstates. Without loss of generality we can 
also assume that 



Bf(K) = B(-K) and w K = -w_ K - 



(69) 



For simplicity we have assumed that we deal with bosonic 
degrees of freedom (note the positive zero-point energy 
in (68)) so that we have 



[B(K),Bt(K' 



JK.K' 



(70) 



The previous commutation relation identifies the opera- 
tors Bt(K) and B(K) as creation and annihilation op- 
erators. Additionally, using the previous formalism, we 
can establish the relations: 

ffB(K) = -cj k B(K) and HB^K) = w K B f (K) (71) 

which shows that the previous creation and annihilation 
operators are eigenoperators of the superoperator H with 
the corresponding eigenvalues ±ojk- In equilibrium all 
operators are related to the Langevin force by a linear 
transformation, and thus, we can express |f) as an ex- 
pansion of the form 



|f)=^c| c B s (K)|^( WK )), 



(72) 



K 



where the kets, |^(wk)), are the eigenfunctions for a 
specific partition of the total system "s" and B S (K) 
are its annihilation (creation) operators. For example, 
the previous expression is valid for the eigenfunctions 
and corresponding eigenoperators of the total EM field- 
matter-bath system as well as for each isolated parti- 
tion containing the bath, like the matter-bath system. It 
is worth mentioning that, in this quantum field theory 
approach, the Hilbert spaces where the eigenfunctions 
and the creation and annihilation operators are defined 
are independent: B S (K) acts on the Fock states defining 
the number of excitations in the field and the eigenfunc- 
tions Ivk(wk)) of satisfy Heisenberg equations of motion. 
From the previous properties and the orthogonality of the 
eigenfunctions one has 



cf,B s (K) = (^( WK )|f). 



(73) 



Using the formal expression for the Green's operator ( 67 1 , 



the eigenbases corresponding to the total and the matter- 
bath Hamiltonians ( 72 ) and the relation ( 73 ) we can take 
the inner product of G??(u>) with |f ) on the right to arrive 
at 
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G N (u) |f ) = [H T 



H 1 



<'+^i^(^))(^(^KOif: 



L K 



K' 

|yK(^K))(yK(^ K )l 

w — uj^ + irj 



L K' 



(74) 



which shows that, for a linear system, the Green opera- 
tor is defined only on the Hilbert space of the eigenfunc- 
tions (i.e. it is not intrinsically quantum) and its poles 
are the mode frequencies of total EM field-matter-bath 
system while its zeros are the mode frequencies of the 
matter-bath system. This result is in agreement with 
the conclusion of Ford et al. [5] . We could have reached 



a similar conclusion by looking directly at Eq. ( 55 ) and 
saying that the mode definition (in this case the modes 
of total system) comes from the only non-trivial solu- 
tion we obtain by imposing P no isc = 0, which implies 



that Det 



G~ N \w) 



(a similar argument can be made 

for the matter-bath system as well) [2]. Note that the 
electromagnetic Green operator is not the only operator 
with which we can retrieve information about the rele- 
vant mode frequencies. Since all parts of the total system 
are coupled, in equilibrium, we could have worked with 
the susceptibility of any one of the other components. 
For example, by working with the matter field we find 



1 + Gjv • xn 



Xn 



IP) = T 



/V 



(75) 



The complex operator T^ 1 in front of the polarization 
field is the the equivalent of the inverse Green operator 
for the EM field. By using (56) T^, 1 can be written as 



1 N 



[1 - G • xn] ■ Xn 1 = G • G^ 1 ■ Xn ' ■ (76) 



From the above equality it is clear that DetT Ar 1 (w) = 
implies 



Det 



Go 



= and/or Det 



G^ 1 



= and/or Det [x^ 1 ] = 
(77) 

where the eigenmodes of the total EM field-matter-bath 
system can be singled out among all possible solutions 
because they are the only ones that depend on the posi- 
tion of the objects. 



C. Sum over modes derivation of Eq. (42 1 



Without taking the continuum limit for the bath os- 
cillator spectrum the operator Gjj 1 ^) is self-adjoint and 
allows for real eigenvalues corresponding with the modes 
of the EM field-matter-bath system. The eigenvalues are 
the solutions of the following characteristic equation: 



Det G)^(w) 



0. 



(78) 



Additional information about the medium can be ob- 
tained from the equation 



Det 







(79) 



which provides the modes in the absence of EM field cou- 
pling. Using the argument principle and a derivation sim- 
ilar to the one of the previous section we can therefore 
say that the energy, no matter what geometry it has, is 
given by 



£(L) 



dz_ 
2ni 



u(z)<9 z lnDet G^{z) 



const. 



(80) 



where, the path C encloses the right part of the complex 
z-planc in the positive sense an infinitesimal distance, S, 
to the right of the imaginary frequency axis. The the ex- 
tra term is due to the energy of the modes of each medium 
and is constant in the sense that it does not change by 
moving one object with respect to the other (no coupling 
without the EM field) . It will disappear when we will take 
the difference with respect to the reference configuration. 

It could be that Det G^ 1 (z) is not a meromorphic func- 
tion because of some branch cuts. As in Schram's deriva- 
tion this problem can be solved by placing the system in 
a "perfect box" and then taking the volume of the box 
to infinity at the end of the calculation (see Fig[IJ). The 
Casimir energy is defined as the difference between the 
energy of a reference configuration (generally, but not 
always, this reference system corresponds with the di- 
electric bodies infinitely separated) and the one we are 
interested in. For outside geometries, the interaction be- 
tween the bodies vanishes at large distances, thus the 
energy at infinity is exclusively determined by the self- 
energies of the dielectric bodies and by the energy of the 
free EM field (see discussion after Eq. (32)). Proceeding 



exactly as in the previous section the finite temperature 
Casimir free energy can be written as 



1=0 



Det 




Det 


G^iib) 



= -fc B T^Trln 



Goo(^) 



(81) 



where we have already taken the dissipative limit and 
used the identity Tr log(M) = log Det(M). The symbol 
Tr indicates the trace on both position and polarization 
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degrees of freedom. We have also assumed that the in- 
tegrand goes to zero sufficiently fast for \z\ oo (high 
frequency transparency) so that the contribution com- 
ing from the integration over the path placed at infinity 
vanishes and also allows us to neglect the surface terms 
after a partial integration is performed. As before an 
equivalent expression in terms of real frequencies can be 
obtained recovering Eq. (42). It is worth noting that 



we would arrive at the same expression for the Casimir 
energy by working with the Green function of the polar- 
ization field. If we write 



F(L) = k B T^\n 



1=0 



Det 




Det 





(82) 



and use Eq. ( 76 ) it is also not difficult to show that 
T 1 — G ■ G 



fc 1 • r 1 



(83) 



so that we recover Eq.(81 1. 



V. SUM OVER COMPLEX "MODES" 

Despite the striking similarity of the final expressions 
for the Casimir energy with and without dissipation they 
are physically distinct from the sum-over-mode stand- 
point, namely, in the dissipative case one must include 
the modes of the bath to obtain the right result. In 
the limit when the spectrum of the bath forms a contin- 
uum, the operator G becomes complex and must fulfill 
the "crossing relation" 



(84) 



where £ is a complex frequency. For this case G is no 
longer Hermitian but, in spite of this complication, the 
Green operator can still be spectrally decomposed as 



6(u) = Q(w) • ]T 



K 



to — cjk 



(85) 



where cjk are now the complex "mode" frequencies (res- 
onances) of the total system and Q{oS) is an entire oper- 
ator (no poles). Passivity requires Imc^K < 0. This is 
just the generalization of (74| ). In order to handle this 
non-self adjoint operator we have introduced the eigen- 
states of the Green eigenvalue equation, I^k), and the 
adjoint-cigcnstatcs (v? K |, they are not simply related by 
Hermitian conjugation [4"B] , i.e. (c^ K | ^ I^k)*- These 
eigenstates are biorthogonal as follows [4*p] 



(Vk'Ivk) = £k',k- 



(86) 



The corresponding eigenvalues, i.e. the frequencies uk, 
can be obtained from the solution to 



Det G-\uj) 



0. 



(87) 



It is worth stressing again that these frequencies are not 
modes and therefore they cannot be summed to obtain 
the Casimir energy |47j . Indeed, a mode represents a 
state of constant energy which can only be defined for a 
closed system where dissipation does not occur. As we 
showed, a formula like equation ([!]) is applicable only 
if we consider the modes of the total system before the 
dissipative limit is taken. 

One may, however, wonder if it is possible to generalize 
equation ([lj to a dissipative system using these complex 
frequencies. To address this point consider the zero tem- 
perature version of the equation given in ( 42 ) written as 



T{L) = E(L) 



h 
'2?r 



/>oo 

/ duj Im 


TrlogG(u;) 
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Only the poles of the Green operator depend on the dis- 
tance L between the objects and therefore the entire op- 
erator Q drops in the previous difference thus we can set 
it 1 in Eq. (85) and in the following analysis. Because of 



the biorthogonality relation ( 86 ) one can show that 
logG 



E I^K)(<^ K l 1 Og (w 
K 



(89) 



After taking the trace (Tt|^k)(^kI = (^kIvk) = 1) and 
performing an integration by parts it is not difficult to 
show 



E(L) 



h 

27T 



duj wlm [d u log D(cj)] c 



(90) 



where D(to) = JIk^ — w k)- Note that Im [d u log D(lo)] 
is the density of states of the EM field. To proceed we will 
decompose D(w) over its zeros. Because of the crossing 
relation ( 84 ) , for every zero wk there is a complementary 
solution — uj^. If we order the frequencies with respect 
to K so that — = cj_k we can write 



d„D(ij) 



E 

I 

E 

K>0 



1 



1 



1 1 



LJ — Wk W + Wjf 



(91) 



where, together with frequency with non zero real part 
w K rI we also allow for the possibility that some of the 
frequencies lay on the negative imaginary frequency axis 
at the position lo^ = — i^^ = i^-K 4 - Physically, the 
pure imaginary frequencies describe a diffusive process 
the EM field undergoes at low frequency when it en- 
ters the medium. This phenonmenon is understand- 
able in terms of eddy (Foucault) currents [31] . For the 
Casimir effect they play a fairly marginal role at zero 
temperature but their effect is rather important at finite 
temperature [Ml [31] . In the last expression we defined 
wk = w k,. , — *^Ki for mathematical convenience, which 
allows all the frequencies to be collected together. We in- 
troduced a prime to stress that the sum over pure imagi- 
nary frequencies has to be taken with a prefactor one half. 
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In Appendix [X] we demonstrate the following indentity 



/Him 



1 



1 



1 L 



W — CJq LU + LUq 



— Im [/(if + 0+)] Re 



Re [/(w )] 

1 i 



(92) 



where "P stands for the principal part, Re [wo] > 0, 
Im [wo] < 0, and /(w) is real along the positive w-axis 
and analytic in left half of the complex plane except the 
imaginary axis. Also, it is assumed that all integrals are 
convergent. Using the previous identity it is possible to 
rewrite Eq. ( 90 1 in the following form 



E(L) = -Re 



K>0 



^2 ( Wk -^Klog 



A 



(93) 



which must be used instead of Eq.Q to obtain the 
Casimir energy if the dissipative limit is taken from the 
outset. This is the generalization of the formula derived 
in Ref. [35] for the plane-plane geometry, which used 
the Lifshitz formula as the starting point, for arbitrary 
geometries. The equation above clearly shows that sim- 
ply summing the complex frequencies obtained from (87) 
does not give the Casimir energy. 



The two terms in Eq.(93 1 have two well-defined origins: 



the first is associated with the energy of the isolated sys- 
tem due to its own vacuum fluctuations (renormalized 
by the presence of the bath) , and the second term comes 
from the coupling with the bath and is associated with 
the bath induced fluctuations. (For finite temperature 
consideration see Ref. [34] ) 

A real cut-off frequency A has been introduced in 



Eq. ( 93 1 for dimensional reason but from the following 



arguments one can easily see that the value A does not 
affect the previous result. Indeed, to get the previous 
formula we have to perform the integral 



t Jo 



Re 



V 



,100 



(94) 



which appears in an intermediate step. It can be per- 
formed noting that for real frequencies D*(uj) = D(—u>) 
and therefore the following sum-rule holds 



1 

2n 



dutuilm [d u log D(w)] = Im 



K>0 



(95) 

This tells us that the term in Eq.(|93|) proportional to 



log A vanishes identically. However, because of the sum- 
rule in Eq. ( 95 1 , the energy contribution of each single 



A. Fluctuation-dissipation theorem vs. the 
sum-over-mode approach 

In the spirit of Casimir's first derivation we have ob- 
tained the Casimir energy using a "field" point of view, 
i.e. summing over the (free) energy of the modes of the 
system. In light of the previous discussion we under- 
stand now that a single mode of the system does not 
refer only to the field subsystem but to the coupled EM 
field-matter-bath system as a whole. In 1955 the ap- 
proach followed by Lifshitz gave a more matter-centric 
view of the Casimir energy regarded as the result of cor- 
relations between the fluctuating currents inside the bod- 
ies. This last approach relies heavily on the application 
of the fluctuation-dissipation theorem (FDT) which can 
be generally stated as a relationship between the sym- 
metric correlation function (quantifying fluctuations) of 
some operator and the imaginary part of the that oper- 
ator's susceptibility tensor (describing dissipation): 



I ({A(w), A(a/)}) = 47r^ImxM<5(w + J) (96) 

where the curly brackets indicates the anti-commutator 
(the symmetric correlation function). This theorem ap- 
plies to both linear and non-linear systems although in 
this last case the calculation of the susceptibility may 
become difficult PHI [SU] ■ The smallness of the external 
perturbations allow us to expand the time evolution of 
the operator up to the first order in perturbation the- 
ory (this expansion becomes simple and exact for linear 
systems). The beauty of the theorem is that it applies 
to a generic operator. A clear example is given by Eck- 
hart in Refs.|441 143] where he calculates the correlation 
function of the EM field coupled to a system of currents 
at thermal equilibrium. The correlation function can be 
either calculated by solving Maxwell's equations sourced 
by random currents and then using the FDT for the cur- 
rents (FDT of the second kind) , or since they are linearly 
coupled to the field, by directly applying the FDT to the 
EM field (FDT of the first kind). 

The "field" and the "matter" point of view seem so dif- 
ferent that they are often described as two distinct meth- 
ods to derive the Lifshitz formula. The sum-over-modes 
approach is generally seen as a technique of limited use 
which does not apply to dissipative systems. In the pre- 
vious paragraphs we demonstrated, however, that this 
depends on the modes we are considering. Here we show 
how the previous considerations connect with the FDT 
and ultimately with the "matter" -point-of- view deriva- 
tion of the Casimir energy. 

First, we would like to clarify a simple but common 
misconception: the susceptibility of a system without 
dissipation is a real function. As an example consider 
the function 



"mode" generally cannot be uniquely defined 



(97) 
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which is the susceptibility for a bulk of matter described 
by dissipation-free oscillators. Before the implementation 
of boundary conditions for the response of this matter, its 
susceptibility (97) is truly real valued and time-reversal 
invariant. However, once boundary conditions have been 
specified for the matter's response to external fields the 
theory of distributions gives the following relation 



where wj£ are the eigenfrequencies of the total EM field- 
matter-bath system. The previous expression clearly re- 
lates to the density of states of the total system and es- 
tablishes a connection between the correlation of the EM 
field as given by the FDT and the modes of the total 
system over which one has to sum to get the Casimir 
effect. 



11TUJ„ 

2wn 



w„) T% + w„)] (98) 



where the sign of the imaginary part mathematically de- 
pends on the way the pole is treated. For example, pick- 



ing both upper signs in ( 98 ) gives retarded boundary con- 
ditions (one could also easily implement time-ordered, 
advanced, and anti-time-ordered response by picking the 
appropriate combination of signs). This is particularly 
clear in a time dependent approach to the system dy- 
namics. The necessity of an imaginary part can be phys- 



ically understood by looking at ( 97 1 as the susceptibility 
of a harmonic oscillator: ujq is the frequency at which 
the oscillator can radiate and / or absorb if set to interact 
with an external field. The fact that in (97) no imagi- 



nary term appears explicitly indicates that the system is 
closed and satisfies time reversal symmetry. 

Let us consider now the application of the FDT to the 
electromagnetic case. The symmetric two-point function 
for the electric field satisfies the following proportionality 
relation in equilibrium 



({E^u, r), E t (u', r')}) cx ImG.^w; r, r'). 



(99) 



The previous expression is the starting point for building 
the stress tensor of the EM field which, in Lifshitz's ap- 
proach, will lead to the Casimir force. The poles of Gi t i 
are complex but from the previous sections we know that 
they can be obtained from the limit where the spectrum 
of the bath oscillators becomes continuous. 

Before taking this limit the previous relation can be 
written as 

{{E N ,i(u, r),E Nji (u;'r')}) cx W^j^u; r, r'), (100) 

where Gj\r,y is the Green function for the EM field in- 
fluenced by the matter and a countable set of bath os- 
cillators. Before the application of boundary conditions 
the Green's function is time-reversal symmetric, meaning 
that G N,i.j is real and even in u). However, it does have 
poles, and by demanding retarded boundary conditions 
we find [3D] 

Im NtiJ (w; r, r ) = Im ^ _' J fc 



K 



^KJ 



ERes^jv,i,j(wK; r, r') r / - ,■ 



(101) 



VI. APPLICATIONS 



Despite the fact that Eq. (81 1 contains all of the nec- 
essary information to calculate the Casimir energy it is 
not the simplest expression one can use for practical pur- 
poses. Similarly, Eq. ( |93[ ), although rich in physical 
understanding is generally inadequate for actual evalu- 
ations. In this section we illustrate some further simpli- 
fications that can be made in some special cases. 



A. Energy of a harmonic oscillator in the 
electromagnetic field. 

Before to approach the calculation of the Casimir en- 
ergy for configurations with multiple bodies it is instruc- 
tive to examine the case when only one body is immersed 
in the EM field 

Due to their mutual interaction, the energy of the EM 
field and the body are shifted from their unperturbed 
values. From a dynamical viewpoint the object's suscep- 
tibility undergoes the process of dressing of by the EM 
field which can be understood as if part of the EM field 
energy forms a cloud of virtual photons around the ob- 
ject [51]. This interaction energy is defined by the total 
energy of the system less the energy of the free EM field 
and the self-energy of the object. From Eq.(81 1 we have 



E 



, f00 dw , 
-n / — coth 

27T 



2k B T 



Im 



Trlog 
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a- 1 ■ Go 



(102) 

where a is the body's bare susceptibility. The Green 
function of the total system is given by 



G 1 = G 1 — a =>■ G — Gq 



1 - a ■ G 



(103) 



and therefore, the previous formula can be written as 



E 







/ — coth 


Jo 2^ 


2k B T 



Im[Trloga d (a;)] (104) 



where the dressed polarizability is defined as 



ad 



I -a- G 



(105) 



This susceptibility describes the dynamics of a body ac- 
counting for the backreaction from the free EM field (see 
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also Eq.(76) and Ref. [52J. Let us consider the sim- 
ple case of a small object with bare dipolar polarizability 
given by 



r a r 



4-> 



Wp{r)p{r') 



(106) 



where p(r) = p(r — r ) is the object's form factor, a 
normalized function (a Dirac-delta in the point-like limit) 
peaked at the position ro of the object. One can show 
that 

(r|d • Go • d|r') = *ct ■ (#„) • ^(r)p(r') (107) 
where 

(#„}(w) = /d 3 r 1 d 8 rap(ra)^o(w,ri ) r2)p(ri) (108) 
and since 

^o(",ri,r 2 ) = J ^^ (k;oj)e^-^ (109) 
we can then write 



(110) 



where p(k) is the spatial Fourier transform of the form 
factor. With similar steps one can prove that 

(rMr') = [1 - # • (^o)]- 1 • Wp(r)p(r') 

= W d p(r)p(r'). (Ill) 

The details of the quantity a <j(w) will depend on the 
dynamics of the polarization field inside the object and 
on the form factor. For an isotropic object we have that 



and in this case Eq.(104) becomes 



E = -3h 
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where the quantity 
V 



J d 3 rp(if 



(112) 



(113) 



(114) 



roughly coincides with the volume of the object. Note 
that even if the bare polarizability does not have any in- 
trinsic dissipation (as it would be for example in the case 
of an atom) the interaction with the free field introduces, 
in addition to a (Lamb) shift of the resonance frequen- 
cies, a radiative damping due to the back reaction of the 
field on the particle. For cases when the object possesses 
intrinsic dissipation (as in the case of a nano-particle) the 
radiative damping augments the intrinsic dissipation. In 
the case of a small object interacting with the free field 
the isotropy of space requires that the radiative damping 



does not depend on its position. This however, is not true 
in the presence of another object where spatial isotropy 
is broken leading to the well-known position dependent 
line shift of the emission spectrum (Casimir-Polder ef- 
fect - see in the following), and to a position dependent 
modification of the decay rate (Purcel effect). One must 
also take care in choosing the form factor: a Dirac-delta 
function leads to some anomalies in the polarizability, 
like acausal response, connected with the longstanding 
discussion of the Abraham-Lorentz equation (see for ex- 
ample [3"3" ] 153} 153] and references therein). 

Before picking a specific form factor we can make some 
general statements. Let us consider the influence of the 
field on a non-dissipative nano-particle described by a 
harmonic oscillator (similar consideration are also valid 
for atoms in the weak coupling limit). First note that 
the equation of motion for the polarizability including 
the backreaction of the field takes the general form 

(115) 

where mo is the oscillator's bare mass, Kq is it's spring 
constant, and o is it's coupling to the field. The backre- 
action term, ( G o(w)), can be simplified by plugging in 
the specific form for the Green's function. For the follow- 
ing discussion it is advantageous to decompose [G q(oj)) 
into its transverse and longitudinal pieces (G*o(lu)) = 

{Gq(lj)) + {Gq(lj)). The transverse part, describing 
propagating waves, is given by 



^(k; W ) =4™ 2 (T- 



<b> kk 



2 / <-» kk 

=4 ™ f 1 ~w 



2lu 



S(k-u) 



(116) 



(117) 



In the second line above we have explicitly introduced 
retarded boundary conditions to enforce causality. The 
two terms in square brackets have distinct physical conse- 
quences: the principal value term will give a mass renor- 
malization and can contribute higher order time deriva- 
tives in the oscillator equation of motion, and the second 
term, containing the delta function, gives the radiation 
damping. The longitudinal component, which gives the 
electrostatic fields from a collection of charges is 



tf£(k; W ) = - 



47rkk 



(118) 



the physical implications of this term will be discussed 
shortly. 

With the explicit formula for the transverse part of the 
Green's function we find 
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4a; 2 
~3T 



(!„) 



(we will assume an isotropic from factor |p(k)| = |p(fc)|). 
The physical interpretation of the longitudinal compo- 
nent is most transparent before the /c-integrals are per- 
formed. The explicit form is given by 



kk' 



(120) 



which is no more than the effective Coulombic spring con- 
stant for a harmonic restoring force (per charge square) 
between two charge distributions described by p(x). 

As a simple example consider a spherically symmetric 
particle with the following form factor 



T|r-r |' 



P(r) 



(2a 2 



IKk)| ; 



(121) 



where the parameter a gives the "radius" of the particle. 
For this specific example we can evaluate the back reac- 
tion ter ms ex plicitly. Using the symmetries of the inte- 



grand in ( 110 ) we ascertain that Cg q(uj)) is proportional 

to the unit dyad (G*o(ui)) = (Go(w)) 4 T > with frequency 
dependence given by 



(GoH) = 



2ur u, 2 a 2 ( „ / iauj \ \ 7T 

^[ 1 + ^ e ^l 1+erf ^JJJ-6^- 

(122) 

The terms in square brackets give rise to a mass renor- 
malization and radiation damping, and the last term, 
coming from the longitudinal part of the Green's func- 
tion, is the shift of the restoring force due to the Coulom- 
bic interaction between the charges constituting the 
dip ole. By inserting our expression for ("G* ( w )) into 
Eq.dllB) 



a.d{ui) = q 2 [—mui 2 + K — iuj"f{uj)\ 
where m is the renormalized mass 



m = too + 
K is the shifted spring constant 
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(124) 



(125) 



K = K 

and 7(w) is the frequency-dependent dissipation constant 



1 (Lo) = — - — e 



1 +erf 



It is worth showing the point particle limit (a 
the polarizability: 

-l 

ad(oj) = q 2 —muj 2 + K — i-q 2 u) 3 

o 

which suffers from acausality. 



(126) 
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B. Many body interaction 

Let us consider now the interaction between two ob- 
jects described by the susceptibilities xi an d Xi (the gen- 
eralization to more than two bodies will become apparent 
from the procedure). For this case the Green function G 
is then given by 



G~ 



G n 



(Xi +X2)- 



(128) 



With only one object, say \2, the Green function G X2 
would be 

G- X l = Go 1 - X 2 => G X2 = Go ■ [1 -X2 ■ Go}- 1 (129) 

and similarly we obtain G Xl for xi- The Green's function 
can always be written as the sum of a free-vacuum part 
plus a scattered part 

G X2 = G a +G X2 => G X2 = G -[1-X2-Go] _1 -X2-G , (130) 

Which allows the total Green's function to be written as 



G X2 ■ [1 - Xid ■ G x 



Gr\ ' G\ 



(131) 



where xid is the field-dressed susceptibility of body 1. In 
the limit that the bodies are infinitely separated Q X2 — > 
indicating that Goo = G X2 ■ Gq 1 ■ G Xl . This is also clear 
from the energetic stand point since this last quantity 
gives the interaction energy of each single body with the 
energy of the free vacuum (which otherwise would be 
over counted). Collecting the previous results Eq.(81) 
becomes 



E{L) = h 



^Trlog[l- ■ * 2 ' G °. ] 

2tt 1 - xi ■ Go 1 - j(2 ■ G 



= h gTrlog[l - xid ■ G Q ■ X2d ■ G ] (132) 

which is the equivalent of the so-called TGTG formula 
discussed in [55H5H] and then re-elaborated for the calcu- 
lation of the electromagnetic Casimir effect within the 
so-called scattering approach [HI I59H6T] . In this ap- 
proach, as one can see in the previous expression, the 
Casimir energy is calculated starting from the knowledge 
of the susceptibility (scattering operator) of each body 
seen as an isolated scatterer interacting with the e.m. 
field (including backreaction) . Note however that in our 
formulation all the previous operators describe dissipa- 
tive systems ab-initio. This means that each body will 
have both a radiative damping, due to the interaction 
with the EM field, and an intrinsic damping, due to the 
coupling with the degrees of freedom of the bath. 



1. Casimir-P older interaction 

Let us consider as a more specific case the Casimir en- 
ergy between a macroscopic body and and small object. 
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In this case it is convenient to begin with 



E = h 



2tt 



THog[l - a d ■ g x ], 



(133) 



where a.d is the dressed polarizability of the small object 
and Q x describes the scattered part of the macroscopic 
body's green function. The expression for the dressed po- 
larizability was derived in the previous section while the 
one for Q x can be easily found in literature for simple ge- 
ometries like a half-space (see for example [52l R)3"| ) . I f 
we perform a series expansion the logarithm in Eq.( 133 1, 
one can show that the first term takes the form 



Tr[d d -£ x ] = / d z v{v\a d -G x \v) = tr[W d ■ <^J], (134) 

where now "tr" traces over the 3x3 tensors and we have 
defined 

tg^){uj,r ) = J d 3 r 1 d 3 r 2 p(r 1 )V x (^ 7 r ll r 2 )p(r2)- 

(135) 

Proceeding in a similar way one can show that 

Tr[(a v -g x ) n ] = J d*r(r\(a d -g x r\r) = tr[(# d ■ ©)"], 

(136) 

and therefore we can re-sum the series to finally get 
E = tt /~ ^trtog[l- r )]. (137) 

The previous is an exact expression for the zero temper- 
ature Casimir-Polder energy that includes multiple re- 
flections. It coincides with the expression given in [53] 
which was derived in a slightly different way for a point- 
like form factor. The usual expansion obtained in second 
order perturbation theory can be recovered by expanding 

the logarithm under the assumption |&d(i£)"\(?x/(*£)l <C 
1 



E 



d£ 
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from where one can get the Casimir-Polder force 
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(139) 



In the case of the exact formula the expression for the 
force looks formally similar but with important differ- 
ences. By differentiating Eq.(133l one has 



F = h 



2tt' 



tr 
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where we have defined 



^Tot(^,ro) = 1 - *ctd(i£) ■ (G x )(i£,,r ) 



1 -(G) (if, r ) 



(141) 



which is the polarizability dressed by the total field (free 
vacuum plus reflected part). This quantity is position 
dependent since its parameter (essentially damping rates 
and resonance frequencies) are determined by the posi- 
tion of the object with respect to the surface. 



VII. CONCLUSIONS 

The role of dissipation in Casimir physics is not yet 
well understood and is at the center of a ten year long 
controversy. This has stimulated several authors to look 
back to the derivation of Lifshitz formula to understand 
if dissipation is properly taken into account. Most ap- 
proaches implement this analysis within the framework of 
the system+bath paradigm which is at the center of the 
open quantum system theory. Up to now Lifshitz's the- 
ory has overcome this careful scrutiny emerging stronger 
than before. 

At first sight Casimir's derivation, based on an ap- 
proach that focused on the modes of the EM field, looks 
rather different from Lifshitz's, and seems incapable of 
being applied to dissipative systems. In this paper we 
have analyzed this point in detail showing that the main 
difficulty is in the concept of mode itself. The complex 
frequencies that we are accustomed to calling "modes" 
in dissipative systems result from a delicate limiting pro- 
cedure. 

We have elucidated this by implementing dissipation 
via the system+bath paradigm. If the bath is com- 
posed of a countable set of oscillators one can show that 
the modes of the system are real; it is the step that 
takes the spectrum of the bath oscillators to a contin- 
uum that leads to complex frequencies. This step is nec- 
essary for the introduction of dissipation since it leads 
to an infinite Poincare time. Knowing this, it is possible 
to calculate the Casimir energy following Casimir's ap- 
proach if we include in our calculation all real modes 
and take the continuum limit only at the end of the 
calculation. We showed that following this approach 
we were able to generalize Schram's calculation [7] to 
obtain the Casimir energy between two parallel plates 
made of a dissipative material. Moreover, the applica- 
tion of the remarkable formula derived by Ford, Lewis, 
and O'Connell 2, 3 generalizes our results for a arbitrary 
geometry. Using the knowledge of the modes of the total 
system one is also able to establish a direct connection 
between Casimir's and Lifshitz's approach through the 
fluctuation-dissipation theorem. 

Despite the difficulties introduced by dissipation we 
also showed that a sum-over-mode analysis of the Casimir 
effect in terms of complex frequencies is still possible. Of 
course this requires a modification of Casimir's original 
expression in order to avoid unphysical complex energies 
[47] , The formula we derived shows, as well as Casimir's 
for non-dissipative systems, a strong analogy with the 
energy of a dissipative quantum harmonic oscillator [261 
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Unlike the Lifshitz formula, where the dissipative and 
non-dissipative cases are related by a minor modification, 
the sum over complex modes formula is rather different 
than the analogous non-dissipative one and clearly shows 
how dissipation enters in the calculation of the Casimir 
energy. In particular it emphasizes that modes which 
have pure imaginary frequency and have no equivalent in 
the dissipationless case must be taken into account for a 
complete description [53] . 

From the previous results it is easy to establish the 
connection with the scattering approach to the evalua- 
tion of the Casimir effect and, in a very simple way, our 
analysis clearly proves its validity in the dissipative case. 
In this last step we showed that in general together with 
any possible intrinsic damping the calculation must in- 
clude a radiative damping due to the interaction of the 
body with the field. It is interesting to note that in agree- 
ment with the scattering approach philosophy, only the 
radiative damping due to the free EM field must be con- 
sidered, which means that each object is treated as an 
isolated scatterer. In all practical calculations the radia- 
tive damping is generally either neglected or phenomeno- 
logically included in the intrinsic damping. 



For simplicity we drop the superscript and the subscript 
and reintroduce them at the end of the calculation. 
The residue theorem gives 
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where an integral around an arc at infinity in the bottom- 
left quarter of the complex plane has been assumed to 
vanish. In the bottom-right quarter we find 
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Because there are no poles in the upper half of the com- 
plex plane we can derive the two integral identities: 
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Appendix A: Sum-Over Poles formula 

Here we give a mathematical demonstration of the 
sum-over complex mode formula presented in the main 
text. Our starting point will be an expression of the form 
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where Re [uq] > 0,lm[cj ] < 0, f(w) is a real function 
along the positive x-axis and analytical in the two left 
quarters of the complex plane except the imaginary axis. 
The quantity ujq may depend on the parameter L. We 
will assume that the initial and intermediate integrals 
are convergent. To ensure convergence a cut-off could be 
adopted which can be relaxed at the end of the calcu- 
lation. The superscript and the subscript on the square 
brackets indicate that we must take the difference of the 
value inside, i.e. 
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and 
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By combining the previous formulas and identities we 
find the relation 



1 



/(w)Im 



1 1 

+ 



+ 
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iujq 



duj = Re [/(w )] 



/ [*(£ - ie) 



(e-ie) a + W g 



(e-^ e )2 + (^r 

If Re [ujq] > we find 



(A7) 
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and when lim^o e |Re [/ — ie)]] j vanishes we can 
write 



7 = Re [/(w )] 



1 



Im [/(i£ + 0+)] Re 



2iujQ 



dt 

(A9) 
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When the imaginary part of the frequency ujq (Im [o>o]=0) 
vanishes we get the expected result 



I = /(wo)- 



(A10) 



As final step both expressions (A9) and ( A12 1 can be 



unified in Eq. ( 92 I noticing that for Re [uq] > 



The case of a frequency ujq with a vanishing real part is 
more subtle. Setting luq = S — i£o (£o > 0) with 6 > 
e in order to keep the pole inside the integration path. 
Sending e — > first we can still use (A9) by replacing 
6 — i^Q. In the limit 



2t(*-i£ ) 



limRe 



and therefore 



/ = Re [/(-i£ + 0^ 



= v 


\ 1 


[e-e \ 


+ )]v 


2£o 1 


e-e \ 



(All) 
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